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The attenuation factor F responsible for the decay out of 
a superdeformed (SD) band is calculated with the help of a 
statistical model. This factor is given by f = (i+rVr s ) -1 . 
Here, Ts is the width for the collective E2 transition within 
the superdeformed band, and is the spreading width which 
describes the mixing between a state in the SD band and the 
normally deformed (ND) states of equal spin. The attenuation 
factor F is independent of the statistical El decay widths Tn 
of the ND states provided that Tn 3> , Ts . This condition 
is generically met. Previously measured values of F are used 
to determine r . 

PACS: 21.16.-n, 21.60.Ev, 21.10.Re, 27.80. +w 

The intensities of E2 gamma transitions within a su- 
perdeformed (SD) rotational band show a remarkable 
feature. The intraband E2 transitions follow the band 
down with practically constant intensity. At some point, 
the transition intensity starts to drop sharply. This phe- 
nomenon is referred to as the decay out of a superde- 
formed rotational band. It is attributed to a mixing of 
the SD states and the normally deformed (ND) states 
with equal spin. The barrier separating the first and sec- 
ond minima depends on and decreases with decreasing 
spin I. Decay out of the SD band sets in at a spin value ±o 
for which penetration through the barrier is competitive 
with the E2 decay within the SD band. The theoretical 
description of this process Jl]-^| uses a statistical model 
for the ND states. The actual physical decay out of the 
superdeformed band is calculated as a function of the pa- 
rameters of the model. Recent experimental data given 
in Refs. ]7|-|Ic| have been analyzed using some of these 
results. 

In the present paper, we take a different approach: We 
calculate the reduction factor F of the intraband tran- 
sition intensity (an observable) directly in terms of the 
spreading width and of the intraband E2 width Tg for 
electromagnetic decay of the SD state. The spreading 
width T* measures the strength of the coupling between 
the SD state and the ND states. We show that F does not 
depend on the electromagnetic decay properties of the 
ND states, provided that the width for statistical El de- 
cay of the ND states is much larger than Tg. Therefore, a 
measurement of F directly yields the spreading width 
and, thus, direct information on properties of the barrier 
separating the first and second minima. This information 
is practically parameter-free and model-independent. It 
is obtained without going through a complete model cal- 
culation of the entire process leading to decay out of a 



SD rotational band. We use data given in Refs. |7| |1C|] 
to deduce values for in 192 ^ 194 Hg and 194 Pb. We test 
our procedure by comparing the values of r§ obtained in 
that same analysis with predictions based on properties 
of the SD band. Our approach and our results differ from 
those of Refs. (2|-||]. We comment on these differences at 
the end of the paper. 

To define the model, we denote the first SD state with 
significant coupling to the ND states during the E2 de- 
cay down the SD band by |0); its energy by E ; the ND 
states having the same spin as the state |0) by \j) with 
j = 1, . . . ,K; their energies by Ej\ and the average spac- 
ing of the ND states by D. We take the limit K — > oo 
in the course of the calculation. The ND states decay by 
statistical El emission. We assume that the total El de- 
cay widths of all ND states have identical values denoted 
by Tn- The coupling matrix elements Voj connect the 
SD and the ND states and are responsible for decay out 
of the SD band. Following Refs. [§-§), we assume that 
the ND states \j) can be modeled as eigenstates of the 
Gaussian Orthogonal Ensemble (GOE) of random matri- 
ces. The spreading width is defined as — 2irv 2 /D, 
where v 2 is the mean square of the matrix elements Voj- 
Typically, we have Tn^s -C D, see Table I. We note 
that in all cases, Ts <Tn. Because of the dominance of 
El over all other transitions, this inequality is expected 
to hold generically. It is fair to expect that decay out of 
the SD band will set in whenever decay time h/Tg and 
mixing time h/T^ of the state |0) are comparable, i.e., 
when is of the same order of magnitude as Ts- Thus, 
we expect that we also have Y* <C D. This is, indeed, 
borne out by the analysis described below. 

The Hamiltonian H of the system is modeled as a ma- 
trix of dimension K + 1 and has the form 

H=( E ° Vaj \ (1) 

To H must be added the diagonal width matrix E with 



E=-(i/2) 







(2) 



The effective Hamiltonian Ti is given by TL = H + S. 

We first investigate an experimental situation which 
would require an energy resolution that is not available 
at present: We study the distribution in energy of the E2 
transition intensity feeding the state |0). Since ^ 0, 
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this state is mixed with the ND states, and the E2 transi- 
tion intensity feeding the state is spread out over a num- 
ber of eigenstates of H. We calculate the ensemble aver- 
age of this intensity distribution. The distribution cannot 
be measured at present because the required resolution 
is of the order of the average spacing D of the ND states. 
However, the calculation defines concepts and yields re- 
sults which are important later on. 

We first consider the case Tn = 0. For T^- <C D, the 
mixing of the SD state with the ND states will cause 
the E2 feeding strength to have a central peak located 
at (or near) Eq, and a number of much smaller and 
well-separated peaks located at (or near) the energies 
Ej . The central peak will obviously be reduced in height 
compared to the case — 0. If observable, this phe- 
nomenon could be referred to as "decay out of the SD 
band" : The peak E2 transition intensity is reduced. The 
reduction occurs even when the ND states are not capa- 
ble of decaying by statistical El emission. If we allow for 
the statistical El decay of the ND states, the picture re- 
mains essentially the same except that all peaks become 
wider. We expect that the reduction of the E2 transition 
peak intensity is not dependent on Tn and is governed 
by the competition between mixing time and decay time, 
i.e., by the ratio F^/Ts- We show presently that this is 
indeed the case. 

To calculate the distribution of transition strength, we 
again consider first the case = 0. The state |0) is fed 
from the next higher SD state. The relative transition 
intensity 1(E) versus energy E is given by 



1(E) = -2 7 [Im (E-E + (i^Ts)- 1 } 7 



(3) 



where Im (or, later, Re) stands for the imaginary (the 
real) part. The symbol 7 denotes the E2 decay ampli- 
tude feeding the state |0), and 1(E) has the expected 
Lorentzian form. Integrating over all E, we obtain 
J dEI(E) = 2tt7 2 . We extend Eq. (|) to the case T l ^ 
by writing it as 



1(E) = -2 7 [Im (E - H)- 1 ]oo 7 



(4) 



where E = E x 1k+i and where Ik is the if -dimensional 
unit matrix. To check the general validity of Eq. (B), we 
transform Ji to diagonal form with left (right) cigenfunc- 
tions ipmi^m) an d complex eigenvalues £ m with m — 
0,...,K. Then, Eq. (g) takes the form 



A" 



1(E) 



-2 7 £ Im [(0|V£) (E - 

7H — 



£ m y l (&)]i, (5) 



which is obviously a sum of Lorentzians. The m th 
Lorentzian is peaked at the energy E = Re £ m . The 
area under the peak is given by 2ir times the probabil- 
ity (0|i/'m)(V'ml > to nnd the state 1°) admixed into the 



in a change of the peak widths of the m resonances and 
the weight factors (0|V , m)(V'm|0)- The total intensity (the 
sum of the contributions of all m resonances) is given by 
27T7 2 , as in the case F^ = 0. It is independent of Tat. 
These results show that we are dealing with a spreading 
width phenomenon, which is not related to the El decay 
of the ND states. 

It remains to average 1(E) over the GOE. We recall 
that the energies Ej are the eigenvalues of a GOE ma- 
trix, and that, correspondingly, the matrix elements Voj 
are uncorrelated Gaussian distributed random variables 
with common mean value zero and common variance v 2 . 
We denote the ensemble average by a bar. Clearly, all res- 
onance structure will be washed out in 1(E), and physical 
intuition and the arguments given above lead us to ex- 
pect that 1(E) has Lorentzian shape and width r$ + IN-. 
This is, indeed, the result of the calculation. To see this, 
we put first E = Eq and write 



[(E - nr\ = -(i/T s ) [1 + S 0Q (E )} 
Here, S o(E) is given by 



(6) 



Soo(E) = 1 - 2ttiJ2woj{D- 1 (E)} j1 w 10 (7) 



j.i 



where Woj — Wjo 



VVWsj V j and 



Dji(E) = (E- Ej + (i/2)r N )^ 7 + inwjowoi 



(8) 



Equations (@J^) show that Soo is an element of a bona 
fide unitary scattering matrix. Aside from the channel 
denoted by zero, there are inelastic channels correspond- 
ing to the emission of El radiation by the ND states and 
leading to the appearance of the term Tn in Eq. Q). We 
have introduced the form Eq. (||) because calculating the 
ensemble average over the diagonal element of a scatter- 
ing matrix defined in terms of the GOE is a standard 
problem in stochastic scattering theory. The ensemble 
average is taken |11 over the distribution of the Voj's 



and Ej's and yields 



S oa (E) = 



E-E a + (i/2)(T s -Tt) 
E-E + (i/2)(Ts+Tl) 



(9) 



For the reader not familiar with the technicalities of 
ref. |pl| , we mention that Eq. ^ can be obtained in 
a simpler way: We replace the ensemble average by the 
substitution E - Ej — * E - Ej + il with I » D. This is 
patterned after Brown's energy averaging procedure p2[ . 
However, we keep E fixed because the GOE average does 
not affect the SD state. The independence of 5oo of Tn is 
caused by the GOE average which replaces the discrete 
spectrum of the ND states by a quasi-continuum. Using 
this result for averaging Eqs. (6) and (4), we find 



eigenstate. A change of the El width T N will result 1(E) = -2j [Im (E - E + (i/2)(T s + T^))" 1 ] 7 . (10) 
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Equation (|10|) shows that the average intensity has 
Lorentzian shape with width Ts + and is independent 
of Tn- The ratio of I(E ) and the intensity for = at 
E = E given by Eq. (3) yields the average peak intensity 
attenuation factor F, 

F = (i + (rVr s ))- 1 . (n) 

The peak attenuation factor F does possess physical 
significance but cannot be measured in practice. How- 
ever, F also determines the reduction of the total in- 
tensity of E2 radiation down from the state |0) and the 
states mixed with it into the next state of the SD band. 
This information is presently available and used below 
to determine F from the data. We have just shown that 
it is not necessary to ask in which way the next lower 
SD state is mixed with the ND states when we ask for 
the total intensity feeding it. It suffices to calculate the 
transition intensity feeding that next lower SD state. 

In the simplest case where ]> = 0, only the inter- 
mediate population of the state |0) is of interest. With 
Tg = 7g, the transition amplitude has the form 

T(E) = 1 (E~E Q + (z/2)r s )- 1 7S . (12) 

We have / dE\T(E)\ 2 = 2tt7 2 : Feeding intensity and 
emitted intensity are, of, course, equal. For T l ^ 0, we 
have 

T(E) = 7 [(E — H) _1 )]oo 7s • (13) 

For T N = 0, the total strength J dE \T(E)\ 2 = 2irj 2 is 
expected to be preserved. This is, indeed, the case and 
can be checked by using the identity 

[(E-W)- 1 )] 00 r s [(E-W*)- 1 )] 00 = 

-2 Im [(E-W)- 1 )]^ ■ (14) 

The identity Eq. (|l4|), in turn, follows immediately when 
we use the form [(E - H)' 1 )]™ = (E - E + (i/2)T s - 
J2j Voj(E — Ej)~ 1 Vjo)~ 1 . Averaging over the ensemble 
does not affect the conservation of total strength. A 
change does occur, however, for Tn ^ 0. Then, the iden- 
tity Eq. (|l4|) does not apply and El decay of the ND 
states weakens the E2 decay intensity. The question is: 
By how much? 

To calculate / dE \T(E)\ 2 for T N ^ 0, we use the form 

[(E-W) _1 ]oo = (l/2)(E-^ + (i/ 2 ) r s) _1 [l + S 0O {E)], 
see Eqs. (6,7), and write Sqo = Sqq + Sqq where the upper 
index fl denotes the fluctuating part. Then, 

^W = 7 2 r s x 

(1/4)[(E - Eo) 2 + (l/4)r|]- l (|l + ^l 2 + • (15) 

The entities on the rhs of Eq. ([l5|) are all known: Soo 

is given in Eq. and | S'qo 1 2 ^ s gi ven explicitly in 

Ref. 0. Integration over \T(E)\ 2 then yields the E2 



intensity down to the next SD state. The actual cal- 
culation would be rather involved. It would involve a 
fourfold numerical integration. Fortunately, the explicit 
calculation is not needed since forr N »r s ,r±, the term 

ISqqI 2 can be neglected. This is seen as follows. As 
in any stochastic reaction problem, the fluctuating part 
of the scattering matrix describes those processes where 
the long-lived intermediate resonances (the ND states) 
undergo statistical decay. For Tn 3> rs,r^, such de- 
cay will overwhelmingly lead to El emission. The term 

\Sqq\ 2 describes the statistical decay of the ND states 
back into the SD state and can, therefore, be neglected. 
We quantify this statement by comparing the sum t^ 
of the transmission coefficients for El decay with the 
transmission coefficient to = 1 — |iSoo| 2 for decay back 
into the state |0). For <c D, we have tn = T^/D, 
while to depends upon energy and is given by to = 
T s T l /[(E- E ) 2 + (1/4) (r s +r^) 2 ]. For ND resonances 
a distance D from Eq we have to ~ TsT^/D 2 <C £n- 
The inequality to <C tn remains valid down to distances 

~ ior s <D. 

Neglecting l-S^I 2 , using in Eq. ( |l5| ) for S o the value 
of Eq. (||), and integrating over all E, we obtain 

J dE\T(E)\ 2 = 2tt7 2 F , (16) 

with F given by Eq. ([Tl|). This is our central result: The 
intensity attenuation factor due to decay out of the SD 
band is given by F. 

Equation ( Jig ) has a simple interpretation. Decay out 
of the SD band is a sequential process with two intrinsic 
time scales, the spreading time %/T* for populating the 
ND states from the SD state |0), and the time h/T-^ for 
El emission. The larger of these two times defines the rel- 
evant overall time scale. This is h/VK Thus, the spread- 
ing width rl signifies the effective partial decay width 
for the total El decay out of the SD band. The branch- 
ing ratio for E2 decay is given by F — Ts/(rs + IN-). 
The branching ratio for the total El decay is correspond- 
ingly given by (1 — F). This can be verified by extending 
Eq. (|l^) to the case Tn =/= 0. These facts make it possible 
to determine directly from the intensity attenuation 
within the SD band. We expect T* to increase strongly 
as we move down the SD band. We point out that in 
contrast to the problem studied in the first half of this 
paper (where F contains the spreading width of the 
state fed by E2 radiation), F now contains the spread- 
ing width rl of the state from which the E2 radiation is 
emitted. 

For the population of the ND states \j) from the SD 
band, the state |0) acts like a doorway state. Because of 
the barrier separating the first and second minima, this 
doorway state has similarity to an isobaric analog reso- 
nance, where isospin conservation has the same function 
as the barrier. The two cases differ in that the width Ts 
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of the doorway state |0) is small compared to the widths 
Fn of the ND states. For isobaric analogue resonances, 
the converse situation holds. A situation similar to the 
one studied above is met in vinydelene, a molecule with 
a shape isomer which is analogous to the SD state ]l3| . 
Here, however, the inequality 3> Tn applies. 

To sum up: Our result, Eq. (|TT|) , is obtained by cal- 
culating Sqo and by showing that for Tn 3> rs,r J -, Sq 
is negligible. We have applied this result to data given 
in Refs. |?] |hJ. The measured quantities are F exp , the 
intensity reduction for E2 decay from the state |0) with 
spin I , and the lifetime r of that state. We equate F exp 
with the ensemble-averaged intensity reduction factor F 
calculated above, F oxp = F. With 1 - F = r l /(T s + T l ) 
and (rg + r^) = Ti/t, we have 

T l = h(l - F cxp )/t. (17) 

The decay width Tg is given by Tg = hF CKp /r. These 
formulae entail an error due to statistical fluctuations of 
-Foxp- The fluctuations are caused by the fact that in any 
given nucleus, the level spectrum is discrete. A very con- 
servative estimate of this error (based on the influence of 
the closest-lying ND states) leads to an uncertainty of 
a factor 2 or 3 for I>. The results are shown in Table 
I. We note that measurable decay out of the SD rota- 
tional band sets in whenever T- 1 > 0.1 Ts or so. The 
values of Tg are consistent with approximately constant 
quadrupole moments in the SD bands. Our results sup- 
port the connection postulated in Refs. Jl],f7| between 
and the barrier separating the first and the second min- 
ima. 

In Refs. I^]] the decay by El (E2) emission is cal- 
culated by multiplying the squares of the projections of 
the eigenfunctions of the Hamiltonian Eq. (1) onto the 
ND (the SD) states with T N (r s , respectively). This 
scheme is expected to apply in the regime r l > r N ,r s . 
Motivated by the later results of Refs. jrHloJ, we have in- 
vestigated the regime Tn 3> I^, Tg. Naturally, decay out 
of the SD band follows different rules in the two regimes. 
Our analysis shows that the regime Tn 3> , Tg actually 
applies to the data summarized in Refs. 0-0. 

In summary, we have shown that data on the attenu- 
ation of the transition strength in a SD rotational band 
yield direct information on the spreading width and, 
thus, on properties of the barrier separating the first and 
second minima. This information is practically model- 
independent and parameter-free. 
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Nucleus 


r N 


D 


r s 


F l 


^Hg(12) 


103 


34 


1.16 


0.18 


192 Hg(10) 


103 


30 


0.54 


5.44 f 


194 Hg(12) 


181 


92 


1.44 


0.97 


194 Hg(10) 


184 


79 


>0.47 


> 8.9 f 


194 Pb(10) 


16 


1699 


0.66 


0.11 


194 Pb( 8) 


17 


1549 


0.28 


0.09 f 



TABLE I. The spreading widths F^ deduced from the data 
reviewed in Ref. H for a number of nuclei. The spin values 
of the decaying states are given in brackets. The units are eV 
for D, and 10 -4 eV for Fn, Fs and r . The results indicated 
with f were calculated with estimated lifetimes in Ref. jjjj. 
The total width r = Ts + F l . 
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